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ABSTRACT 

We present observations of the young Supernova remnant (SNR) 
RXJ1713. 7-3946 with the Fermi Large Area Telescope (LAT). We clearly de- 
tect a source positionally coincident with the SNR. The source is extended with a 
best-fit extension of 0.55° ±0.04° matching the size of the non-thermal X-ray and 
TeV gamma-ray emission from the remnant. The positional coincidence and the 
matching extended emission allows us to identify the LAT source with the super- 
nova remnant RX J1713.7— 3946. The spectrum of the source can be described 
by a very hard power-law with a photon index of T = 1.5 ±0.1 that coincides in 
normalization with the steeper H.E.S.S.-detected gamma-ray spectrum at higher 
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energies. The broadband gamma-ray emission is consistent with a leptonic origin 
as the dominant mechanism for the gamma-ray emission. 

35 Subject headings: gamma-ray: observations; ISM: supernova remnants, ISM:individuals:RX 

36 J1713. 7-3946, acceleration of particles, radiation mechanisms: non-thermal 



Introduction 



38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 

51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 



Gamma-ray observations of shell-type supernova remnants (SNRs) hold great promise 
to help understanding the acceleration of cosmic rays (CRs). These particles - arriving at 
Earth mostly in the form of protons - are thought to be accelerated by a mechanism called 



diffusive shock acceleration (Bell 1978 Blandford & Ostriker 1978 Jones & Ellison 1991 



Malkov & Drury 2001) in the shocks of supernova explosions up to energies around the 
"knee" in the spectrum of cosmic rays (~ 10 15 eV). In particular, X-ray and TeV gamma-ray 



observations of young SNRs such as Cas A (Hwang et al. 2004; Gotthelf et al. 2001 Albert 



et al. 2007; Abdo et al. 2010b), or RX J1713. 7-3946 (Koyama et al. 1997 Uchiyama et al. 



2002 Aharonian et al. 2006, 2007) have confirmed the existence of relativistic particles in 



the shock waves. Young SNRs are preferred targets for seeing particle acceleration at work 
since in these objects the shocks are still strong and actively accelerating particles to the 
highest energies. Gamma-ray instruments have the angular resolution to spatially resolve 
some of the closer SNRs. 

RX J1713.7— 3946 (also known as G347.3— 0.5) is a young "historical" remnant suggested 
to be associated with the appearance of a guest star in the constellation of Scorpius in AD393 
by |Wang et aT| ( p97| ). RX J1713. 7-3946 is located in the Galactic plane (at I = 347.3°, b = 
—0.5°) and was discovered in soft X-rays in 1996 in the ROSAT all-sky survey (Pfeffermann 



& Aschenbach 1996). At a suggested distance of 1 kpc (Koyama et al. 1997; Fukui et al. 



2003 Cassam-Chenai et al. 2004) with angular diameter ~ 65' x 55', the size of the shell 
is ~ 20 pc. Its properties are strikingly dominated by non-thermal activity. Its X-ray 
emission shows a featureless spectrum interpreted to be completely dominated by X-ray 
synchrotron emission from ultra-relativistic electrons ( Koyama et al.|[T997 ; Slane et al.||1999 



Tanakaet aL|2008"] ) . The thermal X -ray emission as well as the radio emission are extremely 



faint (Lazendic et al. 2004). Detailed X-ray observations with Chandra and XMM-Newton 
unveiled a complex structure of filaments and knots in the shell of the SNR - in particular 



in the western part (Uchiyama et al. 2003 Lazendic et al. 2004; Cassam-Chenai et al. 2004 



Acero et al.||2009 ). A recent study with the Suzaku satellite extended the X-ray spectrum to 



40 keV, a measurement that enabled the determination of the parent electron spectrum 



in the energy range where the spectrum cuts off (Tanaka et al. 2008). 
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67 RX J1713.7— 3946 is the first SNR for which TeV gamma-ray emission was clearly de- 
es tected emerging from the shell. H.E.S.S. measurements provided the first-ever resolved 

69 gamma-ray emission at TeV energies. The TeV emission closely matches the non-thermal 

70 X-ray emission as demonstrated by Aharonian et al. (2006). The energy spectrum of 



71 RX J1713.7— 3946 has been measured up to ~ 100 TeV, clearly demonstrating particle 

72 acceleration to beyond these energies in the shell of the SNR. 

73 While the non-thermal X-rays detected in the shells of young SNRs are clearly generated 

74 through synchrotron emission by ultra-relativistic electrons ( Koyama et al.|19~97 ), the picture 

75 of the particle population radiating the gamma rays is not so clear. The main argument 



76 


revolves around two main emission mechanisms ( 


Aharonian et al. 


2006 


Katz & Waxman 


77 


2008 


Berezhko & Volk 


2008 


Porter et al. 


2006; Ellison & Vladimirov 


2008 


Morlino et al. 



2009) 



but so far, conclusive evidence for either possibility is still missing. One scenario 

79 suggests a connection of the TeV gamma-ray emission with accelerated protons (CRs) by 

so invoking the interaction of these protons with interstellar material generating neutral pions 

81 (/T°s) which in turn decay into gamma rays. A second competing channel exists in the 

82 inverse Compton scattering of the photon fields in the surroundings of the SNR by the same 

83 relativistic electrons that generate the synchrotron X-ray emission. This channel naturally 

84 accounts for the close resemblance between the X-ray and the TeV gamma-ray images. 



jSeveral ways have been suggested to distinguish between these two scenarios (see e.g. Morlino 
et al.|2009[ ) but one of the most promising seems to be the broadband modeling of the spectral 



87 energy distribution (SED). In this regard, observations of young SNRs with the LAT on 

ss board the Fermi Gamma-Ray Space Telescope (Fermi) are of particular importance since 

89 the LAT detects gamma rays in the energy range that bridges sensitive measurements with 

90 X-ray satellites such as Chandra and XMM-Newton and TeV gamma-ray telescopes such as 

91 H.E.S.S., VERITAS or MAGIC. 



92 2. Observation and Analysis 

93 The Fermi-LAT is a pair-conversion gamma-ray telescope with a precision tracker and 

94 calorimeter, each consisting of a 4 x 4 array of 16 modules, a segmented anti-coincidence 

95 detector (ACD) that covers the tracker array, and a programmable trigger and data ac- 

96 quisition system. The incoming gamma rays produce electron-positron pairs in the tracker 

97 subsystem, which allow a reconstruction of the directions of the primary gamma rays using 

98 the information provided by the 36 layers of silicon strip detectors in the tracker. The energy 

99 of the incoming gamma ray is determined from the energy deposited by the electromagnetic 

100 showers in the segmented Csl calorimeter. The ACD subsystem is used as a veto against 
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101 the great majority of cosmic rays that trigger the LAT. The energy range of the LAT is 20 

102 MeV to > 300 GeV with an angular resolution for events converting in the front part of the 

103 detector of approximately 3.5° at 100 MeV, improving to about 0.1° at 10 GeV (defined as 

104 the 68% containment radius of the LAT point-spread function or PSF). Full details on the 

105 instrument and the on-board and ground data processing are given in (Atwood et al. 2009). 

we The LAT normally operates in a scanning mode (the "sky survey" mode) that covers 

107 the whole sky every two orbits (~3 h). We use data taken in this mode from the commence- 

los ment of scientific operations on 2008 August 4 to 2010 August 4. The data were prepared 

109 and analyzed using the LAT Science Tools package (v9rl6pl), which is available from the 

no Fermi Science Support Center Q Only events satisfying the standard low-background event 

in selection (the so-called "Diffuse" class events) and coming from zenith angles < 105° (to 

n2 greatly reduce the contribution by Earth albedo gamma rays, see Abdo et al. 2009a) were 

in used in the present analysis. We use all gamma rays with energy > 500 MeV within a 

lu 20° x 20° region of interest (ROI) centered at the nominal position of RX J1713.7— 3946 

H5 (a = 258.39°, 5 = 39.76°). We chose a lower bound of 500 MeV for this analysis for two 

lie reasons: Due to the relative hardness of the spectrum of RX J1713.7— 3946 compared to 

H7 the Galactic diffuse background, photons with energies below 500 MeV are not effective in 

us constraining morphology or spectral shape of the source. Additionally, the broadening of 

n9 the PSF at low energies might lead to systematic problems of confusion with neighboring 

120 sources in this densely populated region of the sky. To further reduce the effect of Earth 

121 albedo backgrounds, any time intervals when the Earth was appreciably in the field of view 

122 (specifically, when the center of the field of view was more than 52° from the zenith) as well 

123 as any time intervals when parts of the ROI were observed at zenith angles > 105° were also 

124 excluded from the analysis. The spectral analysis was performed based on the P6v3 version 

125 of post-launch instrument response functions (IRFs) which take into account pile-up and 

126 accidental coincidence effects in the detector subsystems (Rando et al. 2009). The binned 

127 maximum-likelihood mode of gtlike, which is part of the ScienceTools, was used to determine 

128 the intensities and spectral parameters presented in this paper. 



1 ht tp : / / fermi . gsfc . nasa.gov /ssc/ 
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2.1. Background sources 
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We adopt a background model for the region which includes components describing 
the diffuse Galactic and isotropic gamma-ray emission [^J It also includes all point sources 
within our ROI which are identified in the 1FGL catalog ( |Abdo et al. 2010a) except 1FGL 
1711.7— 3944c which is spatially coincident with RX J1713.7— 3946. All 1FGL sources are 
modeled with a power-law spectrum using the flux and spectral index values obtained from 
the catalog. Exceptions are the known pulsars in the ROI which we model with a power-law 
with exponential cutoff spectral model. As the parameters for this spectral model cannot 
be obtained from the 1FGL catalog, we keep the flux, spectral index and cutoff energy of 
the known pulsars as free parameters in the maximum likelihood fits of the ROI. Figure [T] 
shows two maps of the point-source detection significance, evaluated at each point in the map 
(TS map) for the region around RX J1713.7— 3946 using photons with energies > 500 MeV. 
The flux of the source is not permitted to be negative, this is why negative fluctuations are 
not visible. The detection significance is shown in terms of the test statistic (TS) of the 
likelihood fit. The TS value is defined as TS=2(ln Lx/Lq), proportional to the logarithm 
of the likelihood ratio between a point-source hypothesis (Li) and the null hypothesis of 



pure background (L ) (Mattox et al. 1996). The significance contours of the TeV emission 



observed from the SNR by the H.E.S.S. telescope array (Aharonian et al. 2006) are overlaid 



on the maps. Panel (a) shows the TS map characterizing the excess emission found in the 
region around RX J1713.7— 3946 over our background model described above. A significant 
TS value is found within the spatial extent of the SNR but also in several regions outside of 
its shell. 

We identify three regions of excess gamma-ray emission which are likely not associated 
with the SNR but belong to background sources not recognized in the first Fermi catalog 
(1FGL). Due to the longer integration time of our analysis (24 months vs. 11 months in the 
catalog) and the corresponding improved sensitivity, the appearance of additional sources 
in our region of interest is expected. We simply denote these sources with the identifier A, 
B, C. The source positions are shown in Figure [T] and given in Table |2j The locat ion of 
source A is consistent with a weak radio source ( |Lazendic et al.|[20 04). It is further identified 
in an internal update of the Fermi LAT catalog using 24 months of data. Source B is 
only 11' from the catalog source 1FGL J1714.5— 3830c and could be an artifact caused by 
unmodeled emission from 1FGL J1714.5— 3830c if this source were spatially extended as has 



been tentatively suggested by Castro & Slane (2010). 1FGL J1714.5— 3830c is modeled as a 



2 The LAT standard diffuse emission models (gll_iem_v02.fits and isotropic-iem-v02.fits), available at 



http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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point source in the 1FGL catalog. However, the catalog source is spatially coincident with 
the SNR CTB 37A which has an extent in radio of ~ 15' (Green 2004). A detailed study of 
the morphology of this source is in progress but beyond the scope of this publication as the 
exact morphology of the CTB 37A source does not significantly affect the spectral analysis 
of RX J1713.7— 3946. For simplicity we just assume the emission from this region to be 
described by two independent point sources, 1FGL J1714.5— 3830c and source B. The third 
additional background source C shown in Figure[T]may be associated with RX J1713.7— 3946. 
It is very close to RX J1713.7— 3946, located about 35' from the center of the SNR. However, 
it is spatially consistent with a local enhancement of molecular gas, observed via the radio 
emission from the CO (J—l— >0) transition (Dame et al. 2001). Furthermore, we will show 
below (see Table [2]) that in a combined likelihood analysis of the spectra of RX J1713.7— 3946 
and the surrounding background sources the emission from source C is considerably softer 
than the gamma-ray emission from the SNR. In fact, both the spectral index and the intensity 
of the source are consistent with expectations of gamma-ray emission from a small cloud of 
molecular gas. Nevertheless, we cannot reject the possibility that at least part of the emission 
attributed to the additional background source C is originating from the SNR shell. While 
we consider source C an independent point source in our standard background model of the 
ROI, we repeat the spectral analysis with a model without this source and account for the 
difference in our estimation of systematic uncertainties. Panel (b) in Figure [I] shows the 
detection significance map for the region around RX J1713.7— 3946 (E > 500 MeV) with 
our standard background model accounted for. A comparison with the significance contours 
from H.E.S.S. suggests a spatially extended emission from the shell of the SNR rather than 
a single point source. 



185 2.2. Centroid and Angular Extent 

we We study the morphology of the emission associated with RX J1713.7— 3946 with a 

187 series of maximum likelihood fits, comparing the TS value for different hypotheses about the 

las shape and extent of the source. We fitted the extension (and position) of the gamma-ray 

189 emission with a disk of varying radius. The emission is found to be significantly extended; 

wo the best-fit position (RA, Dec = 258.50° ± 0.04° tat , -39.91° ± 0.05° tat ) is consistent with the 

i9i center of the SNR within 0.2° and the best-fit radius is 0.55° ± 0.04°. This size is consistent 



192 with that of the X-ray SNR given in Green (2004) as 1.1° x 0.9° in diameter. To confirm 



193 these fits, we test a single point source at the location of the highest excess in the TS map 

194 within the shell of the SNR. We further test a spatially extended source defined by the 

195 shape of the H.E.S.S. significance contours of RX J1713.7— 3946 and an extended source 
we as a uniform disk of 0.55° radius. Finally, we consider two and three independent point 
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197 sources within the shell of the SNR located at the most prominent peaks in the TS map. 

198 A power-law spectrum with integrated flux (between 1 and 300 GeV) and spectral index as 

199 free parameters is assumed for each of the hypotheses. The detailed setup of the likelihood 

200 fit is identical to the one used for the spectral analysis and described with that analysis 

201 (Section 2.3). Table [I] shows the flux, and spectral index of the tested shape and its TS 



202 value in comparison to the background model. The TS values are suggestive of extended 

203 gamma-ray emission from RX J1713.7— 3946. The H.E.S.S. significance map as well as the 

204 uniform disk have a difference in TS of ATS = 61 or 58 (H.E.S.S. /Disk) relative to a single 

205 point source and a ATS = 43 or 40 (H.E.S.S. /Disk) relative to a set of 3 point sources within 

206 the shell of RX J1713.7— 3946. However, the TS value in a comparison to the background 

207 model for both the H.E.S.S. significance map (TS = 77) and the uniform disk (TS = 79) 

208 are almost identical, demonstrating that we are not sensitive to the detailed shape of the 

209 emission region. For the models of RX J1713.7— 3946 considered, the TS value is expected 

210 to follow a ^-distribution with two degrees of freedom in the case that no source is present 



2n (Mattox et al. 1996) and therefore can be converted to a detection significance of ~ 8.5a for 

212 both the H.E.S.S. template and the uniform disk model. The positional and the angular-size 

213 coincidence with the X-ray and TeV gamma-ray emission strongly favors an identification of 

214 the LAT source with the SNR RX J1713.7— 3946. 

215 Fig. [2]shows a series of LAT gamma-ray counts maps of the sky surrounding RX J1713.7— 3946. 

216 We choose an energy threshold of 3 GeV for these maps, higher than the analysis threshold 

217 of 500 MeV, to enhance their resolution. The counts maps are smoothed with an 0.3° wide 

218 Gaussian kernel. This width corresponds to the size of the LAT PSF at 3 GeV (the 39% 

219 containment radius of a 2-D gaussian), averaged over front and back conversions and over all 

220 incident angles. Locations of 1FGL catalog sources in the region are marked by squares. Our 

221 additional background sources are denoted by circles and labeled. The black lines again dis- 

222 play the contours of the H.E.S.S. significance map of RX J1713.7— 3946. Panel (a) shows all 

223 counts in the region. The emission coinciding with RX J1713.7— 3946 is faint; the counts map 

224 is dominated by the Galactic diffuse emission as well as emission from 1FGL J1714.5— 3830c 

225 and 1FGL J1705.5— 4034c. Panel (b) shows a residual counts map after subtraction of our 

226 background model. On this panel a clear excess within the shell of RX J1713.7— 3946 is vis- 

227 ible. Panel (c) finally shows the residual counts after subtraction of our background model 

228 as well as the emission from RX J1713.7— 3946 (using the H.E.S.S. significance map as the 

229 template for the spatial extension). The residual counts are consistent with the expected 

230 statistical fluctuations, i.e the region around the SNR is well described by our model. 



231 



2.3. Spectral Analysis 



232 We adopt the spatial extension model based on the H.E.S.S. significance map as the 

233 default model for the analysis of the spectrum of RX J1713.7— 3946. As discussed in the 

234 previous section, the LAT is not able to distinguish between the two extended source models 

235 that we tested. Therefore, we compare the obtained spectrum from the default model to the 

236 results derived from a uniform disk source model and include the difference in the systematic 

237 uncertainty of the spectrum. In the first step of the spectral analysis we perform a maximum 

238 likelihood fit of the spectrum of RX J1713.7— 3946 in the energy range between 500 MeV 

239 and 400 GeV using a power-law spectral model with integral flux and spectral index as free 

240 parameters. To accurately account for correlations between close-by sources we also allow 

241 the integral fluxes and spectral indices of the nearby 1FGL and sources A, B, C (< 3° from 

242 the center of the ROI) to be free for the likelihood maximization, as well as the spectral 

243 parameters of identified LAT pulsars, instead of fixing them to the 1FGL catalog values. We 

244 redetermine in our fit the normalization of the Galactic diffuse emission model, the index of 

245 an energy dependent (power-law) multiplicative correction factor to it, and the normalization 

246 of the isotropic component. This accounts for localized variations in the spectrum of the 

247 diffuse emission in the fit which are not considered in the global model. 

248 For the Galactic diffuse emission, we find a normalization factor of 0.93 ± 0.01 in our 

249 region of interest and a spectral correction factor index of 0.019 ± 0.002 (the positive sign 

250 corresponds to a spectrum that is harder than in the model). The normalization factor for 

251 the isotropic component is 1.17 ± 0.05. These factors demonstrate the good agreement of 

252 the local brightness and spectrum of the diffuse gamma-ray emission with the global diffuse 

253 emission model. Table [2] summarizes the source parameters obtained as results from this 

254 fit. The table includes the spectral parameters and the TS values of all fitted sources. The 

255 flux above 1 GeV obtained for RX J1713.7— 3946 with our default background model is 

256 F 1000 = (2.8 ± 0.7)xl0~ 9 cm" 2 s" 1 and the spectral index is T = 1.50 ± 0.11. Figure § 

257 shows the uncertainty band obtained from this fit. 

258 In a second step we perform a maximum likelihood fit of the flux of RX J1713.7— 3946 

259 in 7 independent logarithmically spaced energy bands from 500 MeV to 400 GeV (using the 

260 spectral model and parameters obtained in the previous fit) to obtain a spectral energy dis- 

261 tribution (SED) for the SNR. The resulting SED is displayed in Figure [3] as black error bars. 

262 We require a test statistic value of TS>4 in each band to draw a data point corresponding 

263 to a 2 a detection significance. This criterion is not fulfilled for the lowest two energy bands 

264 500 MeV-1.3 GeV and 1.3 GeV-3.4 GeV and accordingly we show 95% flux upper limits for 

265 these bands. 

266 In a final step we estimate the systematic uncertainty on the obtained spectral parame- 
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267 ters by repeating the maximum likelihood analysis for several variations of our default model. 

268 Specifically, we varied the source shape template, the background sources, and the model 

269 of the Galactic diffuse emission. The spectral analysis was performed: a) with the uniform 

270 disk shape replacing the H.E.S.S. significance map template; b) with the closest background 

271 source C removed from the model (see also discussion above); c) using a preliminary list of 

272 sources from the 2FGL catalog in development within the LAT collaboration; d) replacing 

273 the standard diffuse emission model by a refined model that is currently being evaluated in 

274 the collaboration for source analysis for the 2FGL catalog (refined with 24 months of data 

275 and with finer gas maps); e) replacing the standard diffuse model by a model based on the 

276 GALPROP code used in the Fermi LAT analysis of the isotropic diffuse emission. The 



280 
281 



277 GALPROP model is described in Abdo et al. (2010c). For e), i.e. the GALPROP-based 



278 model, we considered the various components of the diffuse emission model separately for 

279 which we then individually fit the normalizations in our likelihood analysis. The compo- 
nents are gamma rays produced by inverse Compton emission, gamma rays produced by 
interactions of CRs with atomic and ionized interstellar gas and gamma rays produced in 

282 the interactions of CRs with molecular gas. The model component describing the gamma 

283 ray intensity from interactions with molecular gas is further subdivided into seven ranges of 

284 Galactocentric distance to accommodate localized variations of the CR and molecular gas 

285 density along the line of sight which are not accounted for in the model. 

286 The same model of the isotropic component was used for all model variations a)- 

287 e). From the model variations a) - e) we obtain a systematic uncertainty of +0.08/- 

288 0.10 for the spectral index of RX J1713.7— 3946 and a systematic uncertainty of (+0.6/- 

289 0.7) xl0~ 9 cm~ 2 s _1 for the flux above 1 GeV on top of the statistical uncertainty. The 

290 systematic uncertainty of the derived flux and spectral index related to the uncertainty in 

291 the LAT effective area was evaluated separately. The uncertainty of the LAT effective area - 



292 estimated from observations of Vela (Abdo et al. 2009b) and the Earth Albedo (Abdo et al. 



293 |2009a[ ) - ranges from 10% at 500 MeV to 20% at >10 GeV. The impact on the spectral pa- 

294 rameters of RX J1713.7— 3946 is a systematic uncertainty of ±0.05 for the spectral index and 

295 a systematic uncertainty of ±0.4 for the flux above 1 GeV. The gray band in Figure [3] displays 

296 the superposition of all uncertainty bands obtained in our variations of the default model. 

297 Figure [4] depicts the model variation (b) resulting in the softest spectrum together with the 

298 fluxes in individual energy bands (black error bars) derived for model (b) using the same 

299 procedure as for the default model described above. The range of systematic uncertainty is 

300 particularly important to consider for comparisons of the spectrum to pion-decay dominated 



3 GALPROP is a software package for calculating the diffuse Galactic gamma-ray emission based on a 
model of cosmic-ray propagation in the Galaxy. See |http:/ /galprop.stanford.edu/| for details and references 
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301 gamma-ray emission models which are generally expected to be softer than inverse Compton 

302 dominated gamma-ray emission models. 



Discussion 
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327 
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329 
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331 
332 
333 
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The positional coincidence between the extended gamma-ray emission detected by the 
Fermi-LAT at the position of RX J1713.7— 3946 strongly suggests a physical association be- 
tween the GeV gamma-ray emission and this young SNR. In addition, the region of bright- 
est LAT gamma-ray emission coincides with the northwestern part of the SNR. From CO 
( J = 1 — 0) observations Fukui et al. (2003) and Moriguchi et al. (2005) suggested that this 



part of the SNR is undergoing complex interactions between the supernova shock wave and 
a molecular cloud. This part is also the brightest region in non-thermal X-rays and in TeV 
gamma rays. The match between the locations of brightest emission suggests that the GeV 
emission is also generated by the population of relativistic particles responsible for the TeV 
gamma-ray and non-thermal X-ray emission. 

The origin of the TeV gamma-ray emission from RX J 171 3. 7— 3946 has been a matter 
of active debate (see Zirakashvili & Aharonian 2010, and references therein). There are two 



competing processes potentially responsible for the shell-like TeV gamma-ray emission from 
RX J1713.7— 3946: Inverse Compton (IC) scattering on the cosmic microwave background 
by relativistic electrons (leptonic model) and 7r°-decay gamma rays resulting mainly from 
inelastic collisions between relativistic protons and ambient gas nuclei (hadronic model). It is 
generally accepted that diffusive shock acceleration (DSA) operates at supernova shocks pro- 
ducing high-energy protons and electrons. However, injection mechanisms of supra-thermal 
particles are poorly known so that the current theory cannot tell us about the number of 
relativistic protons and electrons produced at shocks. This makes it difficult to reliably 
predict the levels of leptonic and hadronic gamma-rays. 

The lack of thermal X-ray lines provided a stringent constraint on the gamma-ray pro- 
duction mechanisms. The luminosity of hadronic gamma-rays scales as n^Wp, where hh 
denotes the gas density averaged over the emission volume (where accelerated protons are 
assumed to be uniformly distributed), W p = £-Esn is a total energy content of accelerated 
protons, and ~ 10 51 erg is the total kinetic energy released by the SN explosion. The 



lack of thermal X-ray emission in SNR RX J1713.7— 3946 (Slane et al. 1999; Tanaka et al. 



2008) severely restricts the gas density in the SNR to be small. Ellison et al. (2010) have 



performed calculations of thermal X-ray emission from shocked plasma with non-equilibrium 
ionization in the case of uniform ambient density, following a hydrodynamic evolution with 
which non-linear DSA theory is coupled; they found that the shocked gas densities required 
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for consistency with the hadronic model are nn J$ 0.2 cm -3 . It should be noted that, taking 
i^SN = 2 x 10 51 erg, one needs £ ~ 1 (i.e., extremely efficient acceleration) for = 0.1 cm -3 
and d = 1 kpc. The extremely efficient (more efficient than usually assumed) transformation 
of the supernova kinetic energy into accelerated particles may lead to very low shocked gas 



temperature (Drury et al. 2009), which in turn could change the density requirement. 



The measurements of GeV gamma-ray emission obtained with the Fermi-LAT pre- 
sented in this paper provide new, crucial information about the particle population respon- 
sible for the gamma-ray emission. We have measured the gamma-ray spectrum of SNR 
RX J1713.7— 3946 in the energy range from 500 MeV to 400 GeV and found that the spectrum 
can be characterized by a hard power law with photon index T = 1.5 ± 0.1(stat)±0.1(sys), 
smoothly connecting with the steeper TeV gamma-ray spectrum measured with H.E.S.S. 
Note that the measured gamma-ray spectrum of RX J 171 3. 7— 3946 now covers five orders of 
magnitude in energy, unprecedented for SNRs. 

The hard power-law shape in the Fermi-LAT energy range with photon index of V = 
1.5 ± 0.1 qualitatively agrees with the expected IC spectrum of the leptonic model, as il- 
lustrated in both Figures [3] and |4| If the leptonic model explains the gamma-ray spectrum, 
the Fermi-LAT spectrum is emitted by a power-law part of the accelerated electrons, and 
therefore we can deduce the power-law index of electrons from the measured photon index. 
Using T = 1.5 ± 0.1, we obtain s e = 2T — 1 = 2.0 ± 0.2. The energy flux ratio of the 
observed synchrotron X-ray emission and the gamma-ray emission means that the average 
magnetic field is weak, B ~ 10 fiG (Aharonian et al. 2006 Porter et al. 2006 Ellison et al. 



2010). The maximum energy of electrons is then E e 



20-40 TeV as determined from 



the Suzaku X-ray spectrum (Tanaka et al. 2008). The presence of synchrotron X-ray fila- 



ments varying on yearly timescales ( Uchiyama et al.||2007 ), if interpreted as being due to fast 
electron acceleration and synchrotron cooling, requires B ~ 0.1-1 mG, which is difficult to 
reconcile with the weak average field. Alternatively, the X-ray variability may be caused by 
time- variable turbulent mag netic fields ( |Bykov et aL|2008[ ) which require a smaller magnetic 
field strength. The filamentary structures and variability in X-rays should be attributed to 



locally enhanced magnetic fields in the case of the leptonic model (Pohl et al. 2005). 



As shown in Fig. [3j several groups have previously presented calculations of IC gamma- 
ray spectra. Detailed comparisons between the observed total GeV-TeV spectrum and IC 
models show that none of the previous IC models matches exactly with the data. Some 
additional complications would need to be introduced to realize a better description of the 
gamma-ray data. For example, the shape of the total IC spectrum could be modified if we 
add a second population of electrons (or even multiple populations) which has a different 



maximum energy (see Tanaka et al.||2008 Yamazaki et al. 2009). Yet another way of modi- 
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fying the IC spectral shape is by invoking more-intense interstellar radiation fields, though 



372 this would require substantial increase in the photon density (see Tanaka et al. 2008). 



Even in the case of the leptonic model, it is important to constrain the level of 7r°-decay 
emission at GeV energies by allowing for a hybrid (leptonic and hadronic) model of the 
GeV-TeV gamma-ray spectrum. For proton number index s = 2 (assumed to be same as 



376 the electron number index: see e.g. Baring et al. (1999) for a discussion of why relativistic 



electron and ion indices should be very similar in non-linear shocks), the GeV flux upper 
limit at 1 GeV corresponds to W p < 0.3 x 10 51 (nH/0.1 cm" 3 ) -1 erg for d = 1 kpc, where 
tt-h denotes the hydrogen number density of X-ray/gamma-ray emitting gas. Therefore, the 
leptonic model does not necessarily mean the proton content in this SNR is unexpectedly 
small. 

The GeV measurements with Fermi-LAT do not agree with the expected fluxes around 



1 GeV in most hadronic models published so far (e.g., Berezhko & Volk 2010). Given the 
current models of diffusive shock acceleration, we can discard the hadronic origin of the GeV- 
TeV gamma-ray emission. The proton number index s ~ 1.5 inferred by the LAT spectrum 
is as small as the asymptotic index of s = 1.5 predicted by extremely efficient CR acceler- 



ation (Malkov 1999, see also Ellison & Eichler (1984) for early indications of this limiting 



behavior). Unless this asymptotic index is realized in the shock waves of RX J1713.7— 3946, 
the hard Fermi-LAT spectrum cannot be ascribed to the vT°-decay emission. However, such 



proton energy distribution is not observed in the current models of efficient DSA (Ellison 



et al. 2010). 



392 4. Summary 

393 We have measured the GeV gamma-ray emission from RX J1713.7— 3946 with the 

394 Fermi-LAT. The emission is extended and shows a size that matches the TeV-detected 

395 gamma-ray emission from this SNR. The gamma-ray spectrum for the SNR has been mea- 

396 sured over more than 5 orders of magnitude combining Fermi-LAT and H.E.S.S. observa- 

397 tions. The spectral index in the Fermi-LAT band is very hard with a photon index of 1.5 ±0.1 

398 which is well in agreement with emission scenarios in which the dominant source of emission 

399 is Inverse Compton scattering of ambient photon fields of relativistic electrons accelerated in 

400 the shock front. The dominance of leptonic processes in explaining the gamma-ray emission 

401 does not mean that no protons are accelerated in this SNR, but that the ambient density is 

402 too low to produce a significant hadronic gamma-ray signal. RX J1713.7— 3946 is the first 

403 remnant where the combination with H.E.S.S. data yields spectroscopic measurements over 

404 more than 5 decade in energy that, in contrast to many of the other LAT-detected remnants 
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405 strongly suggests a leptonic origin of the gamma-ray emission. 
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Source morphology 



Flux a 



Photon index TS 



R. A. 2000 Dec. 



Point source 

2 point sources 

3 point sources 



Extended source (H.E.S.S.)j^ 
Extended source (uniform disk) ^ 



1.2 ± 0.7 
0.5 ± 0.5 
1.2 ± 0.9 
0.5 ± 0.5 
1.2 ± 0.2 
0.4 ± 0.3 
2.8 ± 0.7 
3.2 ± 0.7 



1.85 ± 0.31 

1.68 ± 0.41 
2.13 ± 0.41 

1.69 ± 0.41 
2.10 ± 0.28 
1.61 ± 0.31 
1.50 ± 0.11 
1.49 ± 0.10 



18 
20 



32 
77 
79 



257.94° 
257.93° 
257.85° 
257.93° 
257.85° 
259.00° 



-39.75° 
-39.61° 
-39.86° 
-39.61° 
-39.86° 
-39.81° 



258.50° -39.91' 



a E>l GeV, in 1CT 9 cm" 2 s" 1 

6 TS value in comparison to a model with no source at the position of RX J1713.7— 3946. 

C H.E.S.S. significance map is used as a template for the intensity of the gamma-ray emission. 

d A uniform disk with 0.55° radius is used as a template for the intensity of the gamma-ray emission. 
The specified coordinates correspond to the center of the disk. These parameters are the best-fit parameters 
when simultaneously fitting the position and the extension. 

Table 1: Results of the morphological analysis of the gamma-ray emission from 
RX J1713.7— 3946. The integral flux between 1 and 300 GeV and the spectral index are 
the free parameters of the fit and are fitted in the energy range 500 MeV to 400 GeV. 
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Source name 



Flux" 



Photon index Exp. cutoff 



TS 



R. A. 2000 Dec. 



1FGL J1705.5-4034c 


2.1 ± 0.7 


2.16 


± 


0.19 




20 






1FGL J1709. 7-4429 


175 ± 6.4 


1.74 


± 


0.03 


4.46 ± 0.23 


50064 






1FGL J1714.5-3830c 


9.8 ± 1.3 


2.47 


± 


0.09 




228 






1FGL J1716.9-3830c 


1.9 ± 1.1 


2.47 


± 


0.34 




14 






1FGL J1717.9-3729c 


4.9 ± 0.7 


2.34 


± 


0.11 




81 






1FGL J1718.2-3825 


8.4 ± 4.3 


1.64 


± 


0.41 


1.72 ± 0.65 


165 






source A 


1.6 ± 0.5 


2.03 


± 


0.17 




28 


258.84° 


-40.46° 


source B 


4.2 ± 1.2 


2.48 


± 


0.16 




43 


258.71° 


-38.70° 


source C 


2.5 ± 0.7 


2.45 


± 


0.22 




21 


257.47° 


-39.75° 


RX J1713. 7-3946 


2.8 ± 0.7 


1.50 


± 


0.11 




77 







ax 



°in GeV 

c Diffcrence in TS value in comparison to a model with no source at the position of the respective source. 

Table 2: Results of the spectral analysis of the gamma-ray emission in the ROI centered at 
RX J1713.7— 3946. The integral flux between 1 and 300 GeV and the spectral index are the 
free parameters of the fit and are fitted in the energy range 500 MeV to 400 GeV. 
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Fig. 1. — Panel (a): Map of the test statistic (TS) for a point source in the region around 
RX J1713.7— 3946 obtained in a maximum likelihood fit accounting for the background 
diffuse emission and 1FGL catalog sources. Only events above 500 MeV have been used in 
this analysis. H.E.S.S. TeV emission contours are shown in white (Aharonian et al. 2007). 
Rectangles indicate the positions of 1FGL sources in our background model, Several TS peaks 
outside the SNR shell are visible. The 3 peaks marked by circles are added as additional 
sources to our background model (see text). Panel (b): Same map as panel (a), but with 
the 3 additional sources now considered in the background model. 
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Fig. 2. — Panel (a): Counts/sq. deg. observed by the Fermi LAT above 3 GeV in the 
region around RX J1713.7— 3946. The map is smoothed with a 0.3°-wide Gaussian kernel 
corresponding to the width of the LAT PSF at 3 GeV. H.E.S.S. TeV emission contours 
are shown in black (Aharonian et al. 2007). Rectangles indicate the positions of 1FGL 
sources. Circles indicate the additional sources considered in our background model. Panel 
(b): Residual counts after the subtraction of the counts attributed to the background model. 
Panel (c): Residual counts after the subtraction of the counts attributed to the background 
model and to RX J1713.7-3946. 
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Fig. 3. — Energy spectrum of RX J1713.7— 3946 in gamma rays. Shown is the Fermi-LAT detected emission 
in combination with the energy spectrum detected by H.E.S.S. (Aharonian et al. 2007). The green region 
shows the uncertainty band obtained from our maximum likelihood fit of the spectrum of RX J1713.7— 3946 
assuming a power-law between 500 MeV and 400 GeV for the default model of the region. The gray region 
depicts the systematic uncertainty of this fit obtained by variation of the background and source models. 
The black error bars correspond to independent fits of the flux of RX J1713.7— 3946 in the respective energy 
bands. Upper limits are set at 95% confidence level. Also shown are curves that cover the range of models 
proposed for this object. These models have been generated to match the TeV emission and pre-date the 
LAT detection. The top panel features predictions assuming that the gamma-ray emission predominately 



originates from the interaction of protons with interstellar gas (brown: Berezhko &: Volk (20081, blue: Ellison 



& Vladimirov (20081, cyan (solid/dashed): Zirakashvili & Aharonian (2010)). The bottom panel features 



models where the bulk of the gamma-ray emission arises from interactions of electrons with the interstellar 



radiation field (leptonic models), (brown: Porter et al. (2006), blue: Ellison & Vladimirov (2008), cyan: 



Zirakashvili fc Aharonian| (|2010|)). See text for a qualitative discussion of these models. 
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Fig. 4. — Same as Figure [3] but featuring the source and background model which resulted 
in the softest spectrum for RX J1713.7— 3946 instead of our default model. 
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